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ABSTRACT 

We investigate large-amplitude baryon acoustic oscillations (BAO's) in off-diagonal entries of 
cosmological power-spectrum covariance matrices. These covariance-matrix BAO's describe 
the increased attenuation of power-spectrum BAO's caused by upward fluctuations in large- 
scale power. We derive an analytic approximation to covariance-matrix entries in the BAO 
regime, and check the analytical predictions using A^-body simulations. These BAO's look 
much stronger than the BAO's in the power spectrum, but seem detectable only at about a one- 
sigma level in gigaparsec-scale galaxy surveys. In estimating cosmological parameters using 
matter or galaxy power spectra, including the covariance-matrix BAO's can have a several- 
percent effect on error-bar widths for some parameters directly related to the BAO's, such as 
the baryon fraction. Also, we find that including the numerous galaxies in small haloes in a 
survey can reduce error bars in these cosmological parameters more than the simple reduction 
in shot noise might suggest. 

Key words: cosmology: theory - large-scale structure of Universe 



1 INTRODUCTION 

Humans are fortunate that there are baryons in the Universe. Not 
only are we made of them, but baryons are responsible for fea- 
tures in the shape of cosmological power spectra and two-point 
correlation functions that are quite valuable to cosmologists. These 
features are called baryon acoustic oscillations (BAO's), and are 
imprints of acoustic oscillations in the gas of the early universe 
jPeeb les & Yu 1970: Sunvaev & Zel'dovic h 1970: Holtz manll989l : 
lEisenstein & Hu .1998. ; Meiksin. White & Peacock .199Si) . Their 
presence in galaxy clustering statistics provides a standard ruler 
to measure the relation between distance and redshift, and the ex- 
pansion of the universe at late times (Blake & Glazebrook 2003; 
ISeo & EisensteinlboOSh . They provide one of the main tools cur- 
rently proposed for studying the effects of dark energy. 

BAO's have been detected i n modem l ow-re dshift galaxy 
surveys, both in the 2dFGRS llCole et al.l |20()5|) and SDSS 
( lEisenstein e7ai]|2005l ; iHutsil 120061 ; IPercival et al.ll2007ij) . These 



detections were made using the (two-point) correlation func- 
tion, or its Fourier dual, the power spectrum. Conveniently, the 
BAO regime is on large-enough scales that non-linear effects 
are mild. Still, for precision cosmology, these mild effects must 
be understood, and a re the topic of much rec ent research (e.g . 



Springel et al. 20051; Ijeong & Koma tsu 200 d; iHuffetalJ 
Schulz & White 2006; Angulo et al.ll 2007l ; ISeo & Eisensteii 



2006 



2007 



Smith. Sco ccimarro & Sheth 20071^). 



Non-linear evolution of the matter power spectrum tends to 
dampen or smear BAO's. For example. lEisenstein. Sep & White! 



found that large-scale bulk flows and cluster formation 
produce motions that smear out the BAO peak in the correlation 
function, but that these motions are confined to relatively small 
scales of ~ 10 Mpc in Lagrangian space. They argued that 
these motions roughly preserve wiggles on the largest scales of 
the power spectrum, but wipe out wiggles on smaller scales. The 
effects we describe in this paper, using Eulerian perturbation the- 
ory, likely arise physically from the same large-scale bulk motions. 
The attenuation of BAO's can also be understood by considering 
an additive mode-coupling power spectrum, which rises on small 
scales as structure develops, along with a function which attenu- 
ates the linear power sp ectrum on non-linear scales. In the halo 
model (HM, reviewed in ICoorav & Sheth 2002) . this small-scale 
contribution is the one-halo (Ih) term, and comes from pairs of 
objects within single haloes. A qualitative ly similar effect occurs 
in renormalized perturb ation theory (RPT; iCrocce & Scoccimarro! 
I2OO7! ; iMcDonalj 120071) . which is less empirical than the HM, 
and seems more accurate thro ugh translinear scales. For example, 
ICrocce & Scoccimarro! ( l2007h show that, more physically than in 
the HM, the mode-coupling power spectrum in RPT goes to zero 
on small scales. 

In this paper, we show that wiggles exist in off-diagonal terms 
of matter and galaxy power spectrum covariance matrices, almost 
entirely out-of-phase with BAO's in the power spectrum. We inter- 
pret these BAO's in the covariance matrix as manifestations of the 
suppression exacted on power-spectrum BAO's by power on large 
scales. Regions of the Universe with upward fluctuations in large- 
scale power have more-suppressed BAO's. 
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Figure 1. Baryon acoustic oscillations in olf -diagonal terms of the power- 
spectrum covariance matrix, calculated using leading-order perturbation 
theory for the trispectrum. In the top panel, each colored curve is a row 
of the covariance matrix, divided by Piin(fcl)^, where fci is the wavenum- 
ber of that row. In the bottom panel, the curves are additionally divided by 
Pun(k)- The dashed curves show the approximation of Eq. (5). The dot- 
ted curve in the bottom panel shows the approximation without the first- 
derivative term. 



The body of the paper is organized as follows. In Section[2l we 
discuss analytic predictions for the BAO's in the covariance matrix 
of matter. In Section[3] we test the analytic predictions against A'^- 
body simulations, and investigate the detectability of the wiggles 
in the covariance matrix. Finally, in Section |4l we investigate the 
effect of these covariance-matrix BAO's on cosmological parame- 
ter estimation, focusing on parameters directly related to BAO's in 
the power spectrum. We perform this analysis for both matter and 
galaxy power spectra, using the HM framework. 



2 ANALYTIC PREDICTIONS FOR 
COVARIANCE-MATRIX WIGGLES 

The covariance of the matter power spectrum in a sur- 
vey of volume V (neglecting survey-shape effects) is the 
sum of a Gaussian term, which depends on the square of 
the non-linear power s pectrum, and a term involving the 



spe 

matter trispectrum (e.g. [Scoccimarro, Zaldarriaga & Huil 19991 : 
iHamilton. Rimes & Scoccimarrdl2006l. HRS). 



y s.i 



(1) 



where K,! is the volume of shell i in Fourier space (proportional to 
fcf for logarithmically spaced bins), and Tij is the parallelogram- 
configuration trispectrum averaged over shells i and j; 



C{ki, kj)i^j 
T{k,, 



ki , kj 



J d ki d kj 

* s,% S.J 



(2) 



Figure [T] shows the dimensionless, normalized power- 
spectrum covariance ^(fci, fc) = C(fci, fc)y/[P(fci)^P(fc)], for 



various fci's. The wiggles in T are much stronger than in the power 
spectrum itself; the ratio of peak to trough can be nearly a factor of 
ten. However, as we discuss in Sect. |3.l| the amount of noise in a 
measurement of T is also large. 

For Fig. [T] we use Pnn for the power spectrum P, and the 
leading-order (third-order) PT trispectrum for Tij. In the upper 
panel, the curves are shown at redshift z = 0\ the curves in the 
lower panel are independent of power-spectrum normalization and 
redshift. This is because the PT trispectrum involves three powers 
of -Piin(fc), evaluated at different fc's, as does the denominator in 
the definition of T. We use a linear power spectrum from CAMeQ 
( iLewis. Challinor & Lasenbvl boOOl) . with the sam e ACDM cos- 
mology as assumed for the Millennium simulation (Spring el et al.l 



12005) : the Hubble constant Ho 



73 



km 
s Mpc 



, ficDM = 0.205, 



Q.b = 0.045, the scalar spectral index Us = 1, and the rms den- 
sity fluctuation in spheres of radius 8 /iMpc~^ is as ~ 0.9. The 
dashed curve is an approximation for fci <JC k, given by 



T{ki 0,k) 



5038 36 PL(k)k 1 P{Uk)k' 
2205 35 Piin(fe) ~^ 5 Piin(fc) 



(3) 



The dotted curve shows this approximation without the first- 
derivative term. See Appendix |A]for a derivation of Eq. l|3j. 

For fci <^ k, the dominant terms for a ACDM power spectrum 
are the second-derivative term, which produces wiggles 180° out- 
of-phase with the wiggles in Piin(fc), and the constant. The first- 
derivative term shifts the wiggles slightly. The dominance of the 
second-derivative term over the first-derivative term implies that 
the main effect of large-scale power is that of suppressing power- 
spectrum BAO's, not moving them. 

By analogy with the non-linear power spectrum relative to the 
linear power spectrum, it might seem that BAO's in the non-linear 
T could be significantly attenuated on small scales compared to 
the PT prediction. However, leading order for the trispectrum is 
3rd-order, not Ist-order as for the power spectrum; thus, the PT 
trispectrum should be valid into the mildly non-linear regime. At 
z — 0, 3rd-order PT works to 1% for the power spectrum for k < 
0.065 h Mpc~^, while the linear po wer spe ctrum fails at the ~ 1% 
level already at k ^ 0.01 /iMpc"^ jMatsu bara 2007). 

On small scales, the PT trispectrum prediction is expected to 
fail. Probably the most plausible model currently on small scales is 
the halo mode l , the m atter trispectrum of which was worked out by 
ICoora v & hJ feOOl', CH). 

Figure[2]compares the HM and PT predictions for T. We con- 
tinue to use Piin in the denominator of the definition T[k\, k) — 
C{ki, k)V/[P{ki)^ P{k)], for easy comparison of the pure PT pre- 
diction. Alternatively, T could be defined with the halo-model mat- 
ter power spectrum in the denominat or. For th e HM prediction, 
we use the same implementation as in lNevrinck & Szapudi (20d3, 
NS); some details of the implementation also appear in Appendix 
|B] below. At 2; = 0, the three-halo (3h) term contributes signifi- 
cantly on large scales. This cannot be the case at arbitrarily small k, 
where PT should hold exactly. This extra contribution in the trispec- 
trum over PT was noted by CH; it is similar to the shot-noise-like 
one-halo term that cont ributes unphysical power on large scales 
of the power spectrum jCoorav & Shethll200j; ISmith et alj|2003l; 



ICrocce & Scoccimarrol2007l) . However, it could still be that the 3h 
contribution is real at some mildly non-linear scale, such as in the 
BAO regime. 



1 See' http : //camb ■ inf o/| 
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Figure 2. A comparison of different terms in the halo-model pre- 
diction for the normalized power- spectrum covariance T{ki,k) = 
C{ki,k)V/[Piin(ki)'^Plin{k)] at z = 0. We use instead of the 
halo-model matter P in the denominator to show clearly how C(fci , fe) dif- 
fers from the perturbation-theory prediction. On large scales, the Ih and 2h 
terms are negligible, but the 3h term gives a probably inaccurate additional 
contribution. On smaller scales (both larger fci and k) than those shown, the 
Ih term comes to dominate C(fci, k). 



3 TESTING AGAINST iV-BODY SIMULATIONS 

To test the analytic predictions, we ran a large suite of A^-body 
simulations to get a high signal-to-noise covariance matrix. In par- 
ticular, we wanted to test our expectation that the added covariance 
from the HM 3h term on large scales is artificial, and also to test 
the detectability of the wiggles, since their amplitude seems much 
larger than the wiggles in the power spectrum itself. 

We ran 100 particle-mesh (PM) dark-matter simulations of 
box size 1024 Mpc, each with 256^ parti cles on a grid 
256 c ells on a side, solved using the code of lOnedin & Huil 
(11998'). The cosmological parameters we used are the same as 
used for the previous PT plots; for the transfer functions, we used 
CAMB. For the initial conditions (IC's) of the simulations, we 
used the 2LPT (2n d-order Lagrangian Perturbation Theory) code 
(IScoccimarrclfl998l) . run lo z = 49. 2LPT IC's have reduced tran- 
sients over IC's generated from the IZel'dovichl jl970l) approxi- 
mation, and are also more accurate fo r higher-order statistics (see 
ICrocce, Pueblas, & Scoccimarroll200^ and references therein). We 
also ran 400 2LPT simulations to « = 0, 100 of which had the same 
random-number seeds as the PM simulations. 

Figure [3] shows average power spectra from the various 
simul ations. We also sh ow the HALOFIT non-linear power spec- 
trum dSmith et al.ll2003h . The average PM power spectrum departs 
clearly from HALOFIT at A; ~ 0.14 hMpc~^, so this is roughly 
the point to which we trust the PM simulations. The second wig- 
gle at fc 0.12 is attenuated in the PM power spectra relative 
to HALOFIT, but HALOFIT was not developed with particular at- 
tention to BAO's, so the PM simulations may still be trustworthy 
there. The average 2LPT power spectrum exhibits a gradual atten- 
uation in power starting at fc ~ 0.05 /i Mpc^^, but we still analyse 
the 2LPT simulations at higher k for their high-signal-to-noise co- 
variance matrix, keeping in mind this caveat. Scoccimarro (private 
communication) speculates that the deficit in power using 2LPT is 
from the large amount of shell-crossing at 2 = 0, and that smoothly 
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Figure 3. Power spectra of simulations used in the paper, averaged over the 
100 or 400 realizations in each ensemble. The pure-2LPT simulations begin 
a gradual attenuation at rather large scales (fc 0.05 h Mpc~^), while the 
particle mesh (PM) simulations seem trustworthy to fc Ri 0.14 h Mpc~^. 

truncating the initial power for k > 0.3 ft Mpc" ^ would improve 
2LPT's performance. 

To measure the power spectra, we used FFT's run on a 256^ 
mesh with cloud-in-cell density assignment. We corrected for the 
voxel window function by extrapolating the factor between power 
spectra measured from one of the simulations with 256'^ and 1024'' 
grids. We subtracted off the Poisson shot noise of 64 h-'-^ Mpc^ 
even though it is negligible over the range of scales we use. 

Figure |4] shows measurements of T{k) from a few ensem- 
bles of simulations. To measure T, we first measured power spec- 
tra in each simulation, using large bins to reduce noise. The bin 
edges were placed by hand to be approximately at the nodes be- 
tween expected wiggles in T. The first bin runs from < fc < 
0.05 ftMpc"^; the lowest k sampled was 27r/1024 h-^ Mpc 
0.006 /iMpc"^ The edge of the first bin, at 0.05 /iMpc"\ was 
chosen to be about where the wiggles cease being prominent in 
Fig.[T] When trying to detect the wiggles in T, there is a trade-off 
in choosing the edge of the first bin; with increasing k, the signal 
increases with the number of modes, but the wiggles dampen. It 
is possible that the wiggles could be more prominent with another 
choice of bin edge. Subsequent bin edges are shown in the plot, 
where horizontal error bars end. For each ensemble, we calculated 
a covariance matrix from the simulations' power spectra, and found 
T{k) by dividing the covariance in each bin by P{ki)^P{k), where 
P is the power spectrum measured from the simulations, averaged 
in each bin and over all simulations. 

The vertical error bars on the points at fc ~ 0.1 hMpc~^ 
were obtained by jackknife resampling. This method differs from 
spatial-jackknife resampling which is often used in estimating er- 
rors in large-scale structure; in spatial-jackknife resampling, errors 
are estimated by excising small spatial regions from simulations. 
In our case, we form jackknife samples by excluding entire simula- 
tions one-at-a-time from the sample. The error bar width on an A'^- 
simulation point is — 1 times the standard deviation of mea- 
surements of T{k = 0.1 /iMpc"^) from A'^ covariance matrices, 
each formed by excluding one simulation at a time from the sample. 

The outer points at fc ~ 0.07 and 0.12 ftMpc"^ appear with 
cones indicating the robustness of the peak at A: 0.1 /iMpc~^. 
We do not show these outer points with their own error bars be- 
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Figure 4. Measurements of the normalized power-spectrum covariance T 
from ensembles of A^-body simulations. T is measured in large bins to 
beat down noise, with edges roughly where nodes in T's wiggles are ex- 
pected. Only the middle points have error bars, which are typical of all 
three points, since different jackknife samples typically move up and down 
together. The cones coming from the middle points indicate the degree of 
robustness in the peak; see text for details. From 100 PM simulations of box 
size 1024 Gpc, we estimate a detection of the peak at a 5-sigma level. 
The simulations evolved entirely using 2LPT have somewhat lower values 
of T; the 100 2LPT simulations and the 100 PM simulations had the same 
initial conditions. With 300 additional 2LPT simulations, the noise in T is 
reduced, and the peak shape is roughly that of the PT prediction. 

cause there is a high degree of correlation among these three points. 
Typically, the curves for different jackknife subsamples just move 
up or down, and the peak-to- trough contrast hardly changes. The 
half-height of the cone at each ki is the standard deviation (times 
VA^- 1) of the differences [f (0.1 /iMpc"^) - f(k^)\ measured 
from the A'^ jackknife measurements. 

The data in Figure |4] do not indicate any extra contribution 
(such as the HM would give) to T above the PT prediction. Both 
the PT and HM predictions are within the errors of the two left- 
most points for the PM ensemble, and these points may be upward 
fluctuations in T for the first 100 realizations. For the 2LPT en- 
semble, adding 300 more realizations to the sample moves the first 
two data points down, and the last point up. So, we recommend 
the simple PT covariance for the matter power spectrum for k < 
0.1 /iMpc"^. However, the PT covariance must fail at some large 
fc, so for analyses extending to larger k than this, we would rec- 
ommend using the halo-model covari ance, excluding the offending 
3h te rm. This is rather arbitrary, but iNevrinck. Szapudi & Rimed 
( I2OO 6. NSR) showed that the Ih and (less so) the 2h terms domi- 
nate the contribution to Fisher information loss in the halo-model 
matter power spectrum on non-linear scales, and thus the Ih and 
2h terms are important to include. We do not exclude the 4h term, 
since it is where the wiggles we are investigating lie. 

Figure |5] shows T measured with greater resolution. With 
these narrow bins, only for the 400-2LPT ensemble was the level 
of noise low enough to view calmly. The error bars are the square 
roots of the diagonal elements of the T covariance matrix, again 
jackknife-estimated. To reduce the noise, we also measured T us- 
ing the weightings method of HRS. This gives 52 estimates of the 
power spectrum per simulation, by passing two differently phased 
long-wavelength sinusoidal window functions through each vertex, 
face, and edge of the simulation box. When excluding a simulation 
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Figure 5. Higher-resolution (more-finely binned) measurements of the nor- 
malized power-spectrum covariance T from 400 2LPT simulations. The red 
curve is T mea sured using 52 weighting s per simulat ion, according to the 
prescription of Hamilt on, Rimes & Scoccimarrol j2006) . The red curve is 
probably raised from the true level of T in the 2LPT simulation because of 
the beat-coupling which the weightings induce. 

from a jackknife sample using weightings, we excluded all of its 52 
power spectra at a time. 

The reduction in noise from the weightings method is sub- 
stantial, and brings out the shape of T. However, we noticed in 
this and other investigations that the covariance estimated in this 
manner was a bit higher than the covariance from the unweighted 
simulations. This is likely a result of beat coupling to large scales 
that HRS discuss. Beat coupling refers to the extra covariance from 
non-trivial geometry (roughly, edge effects) in a power spectrum 
measured with finite bins. The weightings, even though they are 
minimally invasive, smear power among nearby wavenumbers as 
any survey window function does. When bins have a finite width, 
as they almost always do, this smearing contributes to the covari- 
ance a trispectrum term involving the power spectrum at the 'beat' 
wavenumber (the difference between the nearby wavenumbers). 
We did not use the weightings method for the investigations shown 
in Fig.|4] but a method such as the weightings method would prob- 
ably be necessary to detect covariance-matrix wiggles in real sur- 
veys, and in doing this, beat-coupling should be taken into account. 
However, beat-coupling seems approximately to cause a wholesale 
additive shift in T, so its effects could be negligible if all one cares 
about is the position of wiggles in T. Because of beat coupling, we 
view it as a coincidence that T from the 400 weighted 2LPT sim- 
ulations lies right on the PT estimate; its true level for the 2LPT 
simulations is lower, as in Fig.|4] We suspect that the lower level of 
T in the 2LPT ensemble is related to 2LPT's diminishing accuracy 
on small scales, as shown by the small-scale deficit in its power 
spectra, shown in Fig. [3] 

3.1 Detectability of covariance-matrix wiggles 

We can also estimate roughly how big a survey is necessary to de- 
tect BAO's in T from Fig.|4] In 100 PM simulations 1024 h Mpc*^ 
on a side (with a total volume of 107 hT^ Gpc"^), we detected the 
peak at a 5-cr (4.97-a, to be exact) level. This is the probabil- 
ity that ^(O.OT) < f'(O.l), as estimated by dividing the mean 
of [f'(O.l) - ^(O.O?)] by the standard deviation of this quan- 
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Figure 6. Wiggles in off-diagonal elements of the covariance matrix of 
halo-model galaxy power spectra, as shown for the matter power spectrum 
in Fig.|2] A few different Halo Occupation Distributions are shown, labeled 
by the minimum mass for a halo to host a central galaxy. 



tity (the half-height of the cone). Doing the same on the right 
side of the peak gives an estimate of lOcr for the probability that 
T(O.l) > f (0.12). 

We do not know precisely how the error bars in the difference 
between the peak and the troughs shrink with volume, but we can 
roughly estimate this by comparing the results from the ensembles 
of 100 and 400 2LPT simulations. There are two separate effects 
reducing the error bars: the increase in the number of power spec- 
tra used for the covariance matrix, as well as the increase in raw 
volume. The relationship between error-bar width and volume is 
not necessarily a simple power law. Still, it provides an easy esti- 
mate: this 'difference error bar' at fc ~ 0.7 scales as V^''^^, and at 
k « 0.12 it scales as V°'^^ . Taking the higher of these, an esti- 
mate of the volume of a survey that would give a 1-a detection of 
the peak is 107/5i/°-^^ h-^ Gpc^ = 1.1 h'^ Gpc^ Similarly, 2, 
3, and 4-a detections could be done with 8, 25, and 57 h~'^ Gpc^ 
surveys, respectively. If we conservatively assume that error bars 
oc y/V (as one might first guess), then we estimate that 1, 2, 3, 
and 4-a detections could be done with 4, 17, 39, and 68 Gpc^ 
surveys. 

The wiggles in T would almost certainly be measured with 
galaxies, not matter. Figure [6] shows T of galaxies using the HM, 
using a few different halo occupation distributions (HOD's). These 
HOD's will be discussed in more detail in Sect. 14.11 but for now, 
note that the wiggles in T for galaxies inhabiting smaller haloes are 
approximately like those for the matter, and the wiggles dampen 
somewhat with increasing host-halo mass. However, some of this 
dampening is done by the 3h term, which apparently contributes 
artificially to the HM prediction in Fig. |2] So, there is some rea- 
son to doubt the dampening of the wiggles in T. The situation 
is a bit more complicated here than previously because we di- 
vide by Pg, the galaxy power spectrum for each galaxy sample, 
instead of the linear power spectrum Piin, to get T{ki,k) = 
C{ki,k)V/[P(kifP{k)]. 

Besides the differences in T itself that galaxies would intro- 
duce, several factors would complicate detecting T in reality. Red- 
shift space distortions and complicated survey geometries would 
both need to be taken into account. As the discrepancy with the 
weightings method shows (Fig. O, even a minimal window func- 



tion affects measured covariances. However, in that case, it seems 
that the window function roughly adds a simple constant to T, so 
it is possible that more complicated survey geometries would also 
roughly preserve the shapes of the wiggles. 

There are also factors that would tend to increase detectabil- 
ity. First, as occurs exactly in PT, the galaxy T is likely close to in- 
dependent of redshift, removing one of the usual complications of 
measuring clustering statistics with a huge redshift survey. Also, we 
have not included the fact that a larger survey would allow larger- 
scale modes fci to be accessed, where the wiggles in T{k\, k) are 
more pronounced. And, it is always possible that there exists an- 
other method than we have used to detect the wiggles more effi- 
ciently. 

We have not explicitly investigated the issue of whether the 
wiggles in T could help significantly in fixing the baryon acous- 
tic scale, but it is worth further investigation, since every bit helps. 
One might first guess that it would be easier to measure BAO's in 
a three-point statistic than a four-point statistic such as T. How- 
ever, the particular parallelogram configurations of the trispectrum 
which average to make T depend only on the amplitude of each 
Fourier mode. Real numbers are not only easier to deal with numer- 
ically, but we conjecture that statistics independent of phases have 
smaller cosmic variance, since they are independent of phase cor- 
relations which increase for smaller survey volume. For these rea- 
sons, and for the simplicity of the algorithm of the estimator, these 
configurations of the trispectrum are competitive, perhaps even ad- 
vantageous over using the bispectrum for detecting BAO's. Thus it 
would be interesting to confirm our predictions in data, since this 
would further corroborate our picture of gravitational structure for- 
mation, and, in particular, demonstrate that our understanding of 
BAO's extends to multipoint statistics. 



4 EFFECT ON COSMOLOGICAL PARAMETER 
ESTIMATION 

Baryon oscillations have become quite popular as a prospective 
cosmological probe. Here we consider to what degree BAO's in the 
covariance matrix affect these prospects, both for the matter power 
spectrum and galaxy power spectra. We investigate two parame- 
ters that directly involve BAO's: the baryon acoustic scale, and the 
baryon fraction ft, = nt/{Qb + ^cdm). We do not investigate the 
baryon acoustic scale directly; we investigate the logarithm of the 
sound horizon when the baryons are released from the Compton 
drag of photons (see Eisenstein & Hu 1998, EH). This is propor- 
tional to the baryon acoustic scale, i.e. the location of the peak in 

the linear correlation function. 

We use a Fisher matrix formalism dFished 1 19351 : 
iTegmark. Taylor & HeavensI Il997h to estimate error bars on 
these quantities, in the same way as in NS. The cumulative Fisher 
information in parameters a and /3 over a range of bin indices 
z G 7?. is approximated (assuming a Gaussian likelihood function) 
as 



E 



9 In Pi 1 



dlnPj 
' dp 



(4) 



where Cn is the square submatrix of C with both indices rang- 
ing over TZ. For simplicity, in this paper, we just consider single- 
parameter (unmarginalized) half-error bars; the half-error bar in a 
is cr(a) = 1/VPqq. 

As recommended above in Sect. [3] we use the halo-model 
trispectrum, excluding the 3h term, to compute the matter power 
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Figure 7. The effect of using different covaiiance matrices when analyzing 
matter power spectra to constrain the logarithm of the baryon fraction. Each 
curve starts at a fcmax, and shows how error bars tighten as larger scales 
are included in the analysis. The fiducial covariance matrix we recommend 
(the halo model covariance, minus the 3h term) has been used for the solid 
curves. A no-wiggle power spectrum was used to produce the covariance 
matrices used for the dashed curves. Only the Gaussian variance is used for 
the dotted curves. The wiggles in the covariance matrix reduce CT(ln /j,) by 
7% for fcmax = 0.6. 



spectrum covariance matrix. We do, however, show the results 
when the full HM covariance (including the 3h term) is included; 
it does not make a big difference. We do not use simply the PT co- 
variance because we investigate well into the non-linear regime, up 
to almost k — 1 hMpc^^. 

Figure|7]shows what effect Tij terms in the covariance matrix 
have on estimation of the baryon fraction from the matter power 
spectrum. For this and subsequent figures, we assume a fixed vol- 
ume of 1 h^'^ Gpc^. Each curve starts on the right at a fcmax, and 
shows how error bars in ln/5 tighten as the power spectrum is 
measured over an increasing range of fc. Four curves depart from 
each fcmax curve, calculated using different covariance matrices: a 
purely Gaussian covariance matrix, including just the P{k)^ term 
on the diagonal; a 'HM-3h' covariance matrix, our fiducial covari- 
ance matrix, which uses the full HM except for the 3-halo term; 
a full halo-model covariance matrix; and a 'HM-3h, nowig' co- 
variance matrix, which uses a no-wiggle transfer function of EH 
instead of CAMB to obtain the input power spectrum. We normal- 
ize the 'nowig' power spectrum to match the CAMB power spec- 
trum in amplitude at fc = 10~^ /iMpc^^, and also set it equal to 
the CAMB power spectrum for fc > 130 /iMpc^^, the wavenum- 
ber where the two power spectra cross on scales smaller than the 
BAO regime. Fig.|9](mainly about galaxy power spectra, discussed 
in Sect. I4.lt shows the halo-model matter power spectrum, along 
with its no- wiggle counterpart. 

Going well into the non-linear regime (looking at the blue 
curve, with fcmax ~ 0.6 /iMpc~^), the wiggly covariance ma- 
trix gives about 7% tighter error bars, going to small fc, than 
the non-wiggly covariance matrix. The difference is also 7% if 
a full HM covariance matrix is used for each. The purely Gaus- 
sian covariance matrix clearly gives underestimated error bars here. 
This is because of the 'translinear information plateau;' significant 
correlations arise from halo mass-function fluctuations between 



Figure 8. As in Fig.|2] except showing eiTor bars on the baryon acoustic 
scale (actually, the logarithm of the sound-hoiizon scale). Here, the non- 
Gaussian covariance matrix terms are essentially negligible. The wiggles in 
the covariance matrix enlarge a-(ln /ft) by 2% for fcmax = 0.3. 



S iower-spect rum bins in the range 0.2 < fc/(/iMpc ^) < 0.8 
Rimes & H amilton 2005, 2006, NSR, NS). 

However, in an intermediate regime (looking at the green 
curve, with fcmax ~ 0.2 h Mpc~^), the full covariance matrix ac- 
tually gives tighter error bars than if just the Gaussian term is used. 
Here, the wiggly covariance matrix gives tighter error bars by 4% 
than the non-wiggly one. Starting on larger scales (looking at the 
red curve, with fcmax ~ 0.06 /iMpc~^), the covariance matrices 
including Tij terms again give somewhat larger error bars than if 
only the Gaussian term is used. 

The derivative terms d\nP{k) / d\n ft are the same for all 
calculations in Fig.|7] For this, we calculated linear power spectra 
with CAMB using slightly different baryon fractions, at fixed ag,. 
We then put these varied linear power spectra through our halo- 
model code. The rise of the one-halo term on small scales does 
attenuate the BAO's somewhat compared to the linear power spec- 
trum, but not to the degree se en in A^-body s imulati ons. We thus 



used the meth od described by jPercival et alj ( l2007bl P07), based 
on the work of l Eisenstein et alj i 2007h . The wiggles are attenuated 



by Gaussian window in Fourier space of width 1/(10 /i ^ Mpc). 

Figure [8]is the same as Fig. |7] except it shows error-bar tight- 
ening in a quantity more popular in BAO parameter estimation, the 
baryon acoustic scale (we actually show error bars in the logarithm 
of the sound horizon). In this case, the error bars using the wig- 
gly covariance matrix differ nearly indistinguishably from those 
using the non-wiggly ones; going from arbitrarily large scales to 
fcmax = 0.3, there is only a 2% difference in the error bars. In fact, 
excluding all Tij terms makes at most a 6% difference in the error 
bars. The 'Full HM' curves have been left off this plot, since they 
are indistinguishable from the 'HM-3h' curves. 

In Fig. [8] again we use the same derivative terms 
d In P{k) /d In s for each curve. We estimate 9 In P(fc)/9 In s by 
changing the sound horizon by hand in the transfer function code of 
EH. This is not terribly physical, but it does what we want: chang- 
ing s by hand in the EH code results in power spectra which are 
effectively identical, except that the BAO's are offset. As with our 
investigation of the baryon fraction, we attenuate the baryon wig- 
gles with a Gaussian window. 

It is interesting that the wiggles in the covariance matrix gen- 
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erally result in somewhat narrower constraints in In /;,. Intuitively, 
it makes sense that knowledge of how BAO's are attenuated by 
fluctuations in large-scale power would help to constrain /t, which 
directly affects the BAO amplitude. 

The effect of T is smaller for sound-horizon determination. 
Numerically, this is for two reasons: one, the derivative term 
din P{k)/ da is effectively zero on large scales. Also, the wig- 
gles in the derivative term are well out-of-phase with wiggles in the 
Fisher matrix, C~^. The difference might be greater if the BAO's in 
T were not so close to entirely out-of-phase with the BAO's in the 
power spectrum; that is, T might have a greater effect if large-scale 
power appreciably moved the BAO's. 

4.1 Wiggles in galaxy power spectra 

In the previous section, we showed that BAO's in the covariance 
matrix of the matter power spectrum have a several-percent effect 
on the inferred error bars on the baryon fraction /(,, and a smaller 
effect on error bars in the sound horizon s. Here we investigate to 
what degree these effects propagate through to galaxy power spec- 
tra. 

We use a halo o ccupation distribution (HOD, e.g. 
iBerlind & Weinberg |2002|) to model the galaxy power spec- 
trum and covariance matrix. We build off of the matter trispectrum 
as worked out by CH , and use the satellite HOD as introduced by 
iKravtsov et alj ( I2004 . K04), with its simple Poisson satellite HOD 
moments. See Appendix iBl for details of our implementation. 

In the K04 model, a halo has a central galaxy if and only if the 
halo has a dark-matter mass m > i\fmin. The other two parameters 
in the HOD are Mi, the mass at which, on average, a halo has one 
satellite galaxy, and 7, the slope in the formula giving the mean 
number of satellites in a halo as a function of mass; 

{Ns\m) = {m/Ahy. (5) 

Three parameters give a huge potential parameter space to ex- 
plore. To narrow it, we vary Afmin, and fix Mi /Mmin ~ 30 (which 
K04 suggested as a fiducial value). Observed galaxy power spectra 
at 2; = are generally close to power laws, so to use the most 
physically plausible HOD's, we set 7 at each Afmin by maximizing 
the galaxy power spectrum's straightness. More precisely, we min- 
imized the sum of the squares of the second derivative of Pg (k) in 
log-log space over a range 10~°'* < k/{hMpc~^) < 100. The 
three halo masses 10^^ '', lO^^'^, and IO^^'^Mq have best-fitting 
values 7 = 0.97, 1.23 and 1.68. 

To account for galaxy shot noise, we use the shot-noise-added 
power spectrum (Pg(fc) + 1/rig) for t he Gaussian te rms on the 
diagonal of the covariance matrix (e.g. ICooravll2004h . Here, Tig 
is the number density of galaxies. This treatment of shot noise 
is equivalent to using an 'effecti ve volume' (k) = V/{1 + 
l/[ngPg(fc)]}^ at each k (e.g. iTegmark et al] Il997h . The shot 
noises for the three HOD's used, in increasing order of Mmin, are 
1/n = 6200, 540, and 60 /iMpc"^. For the smallest-Mmin sam- 
ple, the shot noise is negligible up to = 1 ft Mpc~^; for the two 
larger ones, the shot noise equals Pg at fc ~ 0.9 and 0.2 h Mpc~^ . 

Figure |9] shows the power spectra for the HOD's we inves- 
tigate below, along with the matter power spectrum investigated 
in Sect. |4] the CAMB linear power spectrum, and the HALOFIT 
non-linear power spectrum. We also show each power spectrum 
as produced from the EH no-wiggle power spectrum. Each power 
spectrum is divided by the EH no-wiggle linear power spectrum 
for clarity. The wiggles in the matter and galaxy power spectra 
here are only attenuated by the rising, unwiggly Ih term. The only 
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Fi gure 9. Various halo-m odel power spectra used in the paper, divided by 
an feisenstein feHi] i 19981) no-wiggle power spectrum. The colored curves 
are halo-model galaxy power spectra, labeled by the minimum mass of a 
halo which hosts a central galaxy. 



ICP 



b 



10' 





— Matter, HM-3h 

— Galaxies, full HM 

— - Galaxies, full HM, nowig 

Galaxies, no Tij 












































k [/i/Mpc] 


10' 



Figure 10. At fixed volume, en'or-bai' widths on the logarithm of the baryon 
fraction from analyzing power spectra of a few different galaxy samples. 
This is as in Fig. [7] which shows error-bai's for matter power spectra, except 
here, femax is fixed at 0.3 hMpc~'. The wiggles in the covariance matrix 
for galaxies have about the same influence on the error-bars on In as for 
matter, if not a bit more. Also, a significant reduction in error-bar width 
occurs for low A/niin (the minimum halo mass which can support a central 
galaxy). 



place where we implement the more-accurate method of P07 is in 
the derivatives of the power spectrum with respect to parameters, 
which is where it matters significantly for parameter estimation. We 
use the same P07 BAO attenuation for all galaxy and matter power 
spectra, since the attenuation in this model is a function of a single 
parameter, the width of a Gaussian window. It is unclear how to 
change this parameter for different galaxy populations. 

Figures [TO] and [TT] show error-bar tightening in the baryon 
fraction and the sound horizon, in a similar manner as in Figs. |7] 
and[8l Here, we fix fcmax = 0.3 /iMpc~^, but show results for 
three different Afmin's. We include the 3h term in the galaxy co- 
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Figure 11. At fixed volume, en'or-bai' widths on the logarithm of the sound 
horizon s from analyzing power spectra of a few different galaxy samples. 
The wiggly covariance matrix here has a bit greater effect for galaxies than 
for matter; for comparison, see Fig. [8] but notice that the i/-axis is inflated 
here relative to that figure. 



variance matrix since we do not have theoretical arguments and 
measurements (as we do for matter) telling us to discard it. 

At our fixed volume V — 1 h^'' Gpc'', constraints are tighter 
for galaxies with a lower Mmin - There are two causes for this. First, 
galaxies with high A/min are sparser and therefore have more shot 
noise, an effect which wins out over their increase in bias. The shot- 
noise effect can roughly be judged by looking at the dotted curves 
at different A/min's. The difference between the brightest and inter- 
mediate HOD's is mainly from this. 

Second, the dominance of lower-halo (Ih, 2h, and 3h) terms 
grows with Afmin in the HM covariance matrix. This is what causes 
the drastic difference between the middle and faintest samples. 
It also increases the impact of covariance-matrix wiggles for the 
faintest galaxies; for this sample, the difference in error-bar width 
between the wiggly and non-wiggly covariance matrices reaches 
13% for In ft, and 4% for In s. The lower-halo terms attenuate the 
wiggles in T increasingly for brighter galaxy samples, as shown in 
Fig.H 

As Mmin increases, not only are the effects of lower-halo 
terms greater on large scales, but the scale increases (i.e. the crit- 
ical k reduces) where lower-halo terms come to entirely dominate 
the covariance. Thus, if A/min is small, the wiggles in the covari- 
ance matrix can be accessed over a larger range of scales. Roughly, 
if the Ih term comes in on smaller scales in the power spectrum 
(as it does for the fainter galaxy samples), then more cosmological 
information can be measured from the 2h (quasi-linear) term. 

This is similar to what we found in NS; cosmological informa- 
tion is more pristine in regions of the universe in earlier stages of 
structure formation. It seems likely that, again as in NS, imposing 
an Afmax, i.e. excluding galaxies in large clusters, could help sig- 
nificantly to tap smaller-scale information. However, just as in that 
study, we find that our model breaks down when the lowest-mass 
haloes dominate the power spectrum and its covariance matrix. For 
Mn.in < lO"Af0 , using values of 7 fit in the same way as for the 
samples we display, we sometimes get non-positive-definite covari- 
ance matrices. (However, it is possible with slightly different 7's to 
get well-behaved covariance matrices.) 

We attribute this breakdown for small Mmin to our lack 



of accurate knowledge about the power spectrum and trispec- 
trum of haloes (Phh and rhhhh) in the translinear regime. 
ISmith. Scoccimarro & Shettil ( l2007al) showed that Phh is difficult 
to model; Thhhh is likely harder. Our assumption that leading-order 
PT describes the polyspectra of haloes is not unreasonable, but the 
leading order is different for each polyspectrum. Using different 
orders of PT for the power spectrum and trispectrum together is 
inconsistent in the translinear regime, where higher-order correc- 
tions to the power spectrum are significant. For the matter power 
spectrum, and for high A/min, the Ih terms of both the power spec- 
trum, and its covariance, dominate terms involving Phh and Thhhh 
on translinear scales. In these cases, accurate modelling of Phh 
and Thhhh would make a negligible difference in our calculations. 
When the Ih terms are reduced relative to other terms, though, the 
raw Phh and Thhhh are exposed in a regime where they are not 
known accurately. 

Because the Afmin ~ W^^'^ Mq sample is close to the critical 
Afmin ~ 10^^ A^Q where covariance matrices can be non-positive- 
definite for fiducial 7's and A-Zi's, we view the drastic diminution 
in the error bars using this sample with some caution. Still, we con- 
fidently assert that in the HM, reducing Mmin can produce signif- 
icantly smaller error bars than the simple reduction in shot noise 
would suggest. 

Even though there seem to be significant gains in cosmologi- 
cal information when using a fainter galaxy sample, it is important 
to remember that we are holding the volume of the survey fixed. 
We have not considered some important realistic effects, an obvi- 
ous one being that it takes less telescope time to measure bright 
than faint galaxies. 

The question of how to choose a galaxy sample to maximize 
cosmological information is a very important and interesting one, 
but it lies outside the scope of this paper. The answer to the main 
question of this section is this: wiggles in galaxy power-spectrum 
covariance matrices do seem to affect constraints on In ft and In s 
to about the same (if not somewhat larger) degree as for matter. 



5 CONCLUSIONS 

Our main points are the following: 

• In off-diagonal entries in power spectrum covariance matri- 
ces, BAO's exist which appear much stronger than the BAO's in 
the power spectrum. These wiggles are a manifestation of the sup- 
pression which large-scale power does to BAO's in the power spec- 
trum, and originate in the perturbation-theory trispectrum. We give 
a simple analytic approximation to these wiggles in terms of the 
linear power spectrum and its first two derivatives, and check the 
analytical predictions using A^'-body simulations. 

• These wiggles are potentially detectable in current and up- 
coming surveys, but because of the large noise in a covariance ma- 
trix measurement, they are only detectable at a small significance 
level. We estimate that a one-sigma detection could be done with a 
survey of size a couple of /i"'' Gpc^, and that a three-sigma detec- 
tion could require a survey of volume ~ 30 /i"'' Gpc^. 

• The wiggles make a modest difference in cosmological pa- 
rameter error bars from analysing galaxy and matter power spectra. 
For example, using the true, wiggly covariance matrix in estimating 
the baryon fraction results in error bars several percent tighter than 
a no-wiggle covariance matrix. Doing so in estimating the baryon 
acoustic scale has a smaller effect. 

• At fixed volume, including galaxies in smaller-mass haloes 
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provides tighter error bars on parameters such as the baryon frac- 
tion and the baryon acoustic scale than analysing galaxies in only 
large haloes. In the context of the HM, this effect goes beyond the 
simple gains from analysing a sample with smaller shot noise. We 
attribute this to the one-halo term of the galaxy power spectrum 
becoming dominant at smaller scales for small haloes than large 
ones. 

Most of the calculations in this paper made use of our package 
of Python code for cosmology, called CosmoPy. It can be down- 
loaded from http : / /www . if a . hawaii . edu/ cosmopy/ 
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APPENDIX A: WIGGLES IN THE PERTURBATION-THEORY TRISPECTRUM 

The expression for the 3rd-order PT trispectrum includes 12 terms involving the F2 symmetrized kernel, and 4 terms involving the F3 
symmetrized kernel. In this Appendix, all appearances of F2 and ^3 are the symmetriz ed kernels, often de note d in the literature with a 
superscript Expressions for these mode-coupling kernels can be found in, for example, lGoroffetal.l ( ll986l) and lBemardeau et~ak I( l2002h . 
The F2 kernel is 

P2(fei , fea) = i [1 + M + + ^) c + (1 - , (Al) 

where c = cos0i2. In the following, we take the Einstein-de Sitter value of fi = 3/7, even though it has a mild dependence on Qrn', 
fi K, ^Q,^^^^ teouchet et aljl993) . This quantity ^ also appears in the F3 kernel. For calculations, we do use this expression for /i; for our 
assumed cosmology, /i ~ 1.045 x |. The PT trispectrum is 

T(fcl,fe2,fe3,fe4) = 4[i^2(fel3,-fc3)-F2(fe24,-fc4)P3P4Pl3 +CyC.] +6[F3(fc2,fc3,fc4)P2P3P4 +CyC.]. (A2) 

Here, fci3 = fci + ^3 and P13 = P(fci3); in this Appendix, P denotes Pun, the linear power spectrum. 

The abbreviation 'eye' describing the 12 F2 terms could do with some further explanation. The P2 terms come from terms such as 
(5''^^(fei)(5'^-'(fe2)(5'^''(fe3)(5''^'(fe4))c, where the superscript on 5 is the order of perturbation. There are six choices of where to put the '^'^^ 
labels. Each of these gives two P2 terms; 

{<5('>(fcl)5<"'(fe2)<5'''(fe3)'5*''(fc4))c = P2(fel3,-fc3)P2(fe24,-fc4)P3P4Pl3 +P2(fel4,-fe4)P2(fe23,-fe3)P3P4Pl4. (A3) 

The PT trispectrum contribut ion to the power- spectrum covariance can be simplified to the following, given by 
IScoccimarro. Zaldarriaga & HdHl999l) (which we rearrange): 



: / / {I2P1P2 [F3(fel,-fel,fe2)Pl +P3(fe2,-fe2,fel)P2] +8P1-2 [P2(fe2-l,fci)Pi +F2(fel-2,fe2)P2]^} 
J ki Jki 



j£ki_j£k2_ 

(A4) 



Here, fci_2 = fci — fe, and P1-2 — P(|fci_2|). 

A possibly useful intermediate result for the P2 part is, again where c = cos 6112, 

|"(3 - Wc')ki + Tcfe] fciPi + (fci ^ k2) 

P2(fe2-4,fcr)P4 +P2(fc4„2,fe2)P2 = ^ (,2 ^ J2 _ • (A5) 

Where this time c = cos 6ab, the P3 part can be expressed as ( I Valageasll2004t) : 

6^ [a''(10 - 59c2 + 28c'') + 0^6^(10 - 44c^ + 76c*) - 2lb*c^] 
^^('^'-'^'''^^ 126a^ [a^ + b^ + 2a%^il-2c^)] ' ^^^^ 

The main approximation we use is ki ^ kj; we keep terms up to zeroth order in fci (terms in fc2/fci to the powers 2, 1, and 0). 
Substituting P(fe2) — (cfci)P'(fc2) + ^(cfci)^P"(fc2) for P1-2 and averaging over c gives 



^Pfe)-gp'fe)fc, + ±P"(fc,)fc-. 



.2205- 35- ^-^'-^ ■ 5- ^•■^'•■^J ^^^''^ ' ^^^^ 

which is the result in Eq. l|3j- There is no term involving a higher derivative of P(fe) to this order in ki. Another approximation which makes 
a term involving P{ki)P{kj) disappear in Eq. JA7| ( is p^f-)- ^ P{k -) ■ 



APPENDIX B: THE HALO-MODEL GALAXY TRISPECTRUM 

The galaxy trispectrum, excluding the effects of shot noise, has a similar form to the matter trispectrum, which CH worked out. The main 
difference is that matter densities get replaced by galaxy number densities. 

In the HM, galaxies and dark matter are assumed to lie entirely within haloes. Like CH, we assume that the haloes are distributed 
according to leading-order perturbat ion theory (PT), i.e. the linear power spectrum, and the 3rd-order PT trispectrum. This is not quite true 
dSmith. Scoccimarro & Shethll2007al) . but there seems currently to be no better analytic approximation to use. 

The matter power spectrum in the HM is the sum of two terms: 

PmAk) = P"m(fc) + PlUk) = Mt{k, k) + P''-^{k)[Ml{k)f, (Bl) 
where Mf are integrals over the halo mass function n{m). 

M^{ki, . . . ,ka) = I {m/ p)°'u{ki, m) ■ ■ ■ u{ka,'m)b0{m)n{m) dm. (B2) 



Here, p is the mean matter density, bplm) is the /3-order halo bias l iMo, Jing & Whitj|l997l: Scoccimarro et alj|200lL SSHJ), and u{k, m 



is the Fourier- transformed halo matter-density profile, for which we take a Navarro. Frenk & Whitd (Il996h form. In principle, u could be 



complex, but we assume a spherically symmetric form, which forces it to be real. We use a ISheth & Torrnenl i ll 9991) form for the halo mass 
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function. We also integrate over a iBullock etalH200lh concentration distribution at eacii mass, but suppress tliat integral in these equations 
for simplicity. 

Our implementation of the galaxy trispectrum uses the subhalo HOD of K04. There are other parameterizations of the HOD, for example 
a binomial model by SSHJ, but the K04 model gives a simple form for moments of the HOD. The parameters of the K04 model are explained 
in Eq. lO, and the text above it. 

The number of satellites is assumed to obey a Poisson distribution. This gives a simple formula for its factorial moments, which 
iKravtsov et alj ( l2004h found to hold through the third moment. Where A'^i"' is the ath factorial moment of the satellite HOD, 

iVi"'(m) = {NsiNs - l)...{Ns -{a- l))|m) = (iV.|m)" = (m/Mi)"^ (B3) 

The dark matter trispectrum has Ih, 2h, 3h, and 4h terms, for which we refer the reader to CH. When converting these to galaxy 
trispectrum terms, the A/f factors are replaced with G^, defined as 



GS(fci 



) ) 



Ug{k 



i,m 



i— 1 



bf:){m)n{m) dm. 



(B4) 



This is a simplification of expressions bv lSmith. Watts & ShethI | |2007| . e.g. their Eq. 76) for the case of the satellite HOD. Here, fig is the 
mean galaxy number density. For simplicity, we assume that the halo galaxy-number-density profile Ug is real, since we deal exclusively with 
spherically symmetric haloes. Ignoring the second term in brackets and the lower limit of integration, this expression has a nearly identical 
form to A/f in Eq. iB2i . The second term in brackets shows the contribution by multiplets of galaxies with the central galaxy. The central 
galaxy is assumed to lie in the centre of its halo, so in real space, the term including it gets one fewer convolution over the galaxy-density 
profile. 

For example, the galaxy power spectrum is 

^^i^^'-^ -G§(fc,fc)+P""(fc)[G}(fc)]' 



(1)/ 



n(m) dm + P""(fc) 



'gal 



NP{m)ug{k,m) 

figal 



-fei {m)n{m)dm 



(B5) 



We note that expressions for the 2h part of the halo-model galaxy power spectrum in the literature sometimes omit the second term in the 
2h numerator (the number 1, which comes from pairs of central galaxies with each other), even though it dominates the first term in the 2h 
numerator on small scales. However, the Ih term is dominant over the 2h term on small scales, so the central-galaxy contribution to the 2h 
term usually contributes negligibly to the overall Pg. 

For the galaxy power spectrum and trispectrum, we do not change the galaxy-density halo profiles or concentrations from what we used 
for the dark matter. As with the matter covariance matrix, we integrate over concentration para meter at each mass. Galaxy-density profiles 
are almost certainly different from matter-density profiles in reality (e.g. lNagai & Kravtsovll2005i) . but we are unaware of a simple analytic 
alternative to what we used for the dark matter. In any case, these approximations are adequate for our present purposes. 



